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ABSTRACT 

We investigate the impact on low energy flavour mixing phenomena of the renormalisation 

X 

group running of parameters in supersymmetric theories. We have explicitly chosen flavour depen- 
dent supergravity couplings. The renormalisation group equations are displayed in general matrix 
form. We work in a scale dependent basis such that the hierarchy in the quark masses can be more 
easily utilised. Simple analytic solutions are provided in matrix form (for general boundary condi- 
tions) when only the third generation Yukawa couplings are retained in the renormalisation group 
equations. We solve the evolution equations of the KM matrix and the supersymmetry breaking 
soft terms in the squark sector. We suggest a general parametrisation of the soft analytic terms, 
more suitable for non-universal supergravity couplings. We point out that for flavour dependent 
soft terms, there are new contributions to the neutron electric dipole moment . Finally, we consider 
a model with maximal flavour mixing in order to estimate the effects due to the renormalisation 
scale dependence. 
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1. INTRODUCTION 



Twenty years ago, the discovery of the charmed quark established the GIM mechanism and 
naturally explained the approximate flavour conservation in neutral current (FCNC) processes 
[1]. It was then generalised to include the third family with the introduction of the Kobayashi- 
Maskawa matrix. This settled the standard model framework enabling one to calculate, up to 
non-perturbative effects, flavour changing neutral current phenomena as radiative corrections. 
Consequently, the comparison of theoretical expectations and experimental data provides valu- 
able information on possible new physics beyond the standard model. On top of the traditional set 
of parameters which characterise FCNC in the ET-system, more recent studies of £>-meson decays 
have considerably enlarged the FCNC phenomenology. A rich literature [2 — 10] is available about 
FCNC restrictions on supersymmetric extensions of the standard model. Nevertheless, both the 
LEP (and Tevatron) constraints on supersymmetric theories and some fresh insight on sponta- 
neously broken supergravities from superstrings have encouraged a recent revival of this subject 
[11-14]. 

The basic supersymmetry induced FCNC effects are produced by the supersymmetric ana- 
logues of the standard model loop diagrams for neutral current processes, where gauge bosons are 
exchanged. Quarks and vector bosons are then replaced by squarks and gauginos. It is convenient to 
separate such graphs into two kinds: (a) Loops that generate FCNC effects in the standard model 
with the VF-bosons replaced by charginos. (b) Flavour preserving loops in the standard model 
with gluons, photons and Z-bosons replaced by gluinos and neutralinos. If equal-charge quark 
and squark mass matrices are not diagonal in the same basis, their couplings to neutral gauginos 
can induce FCNC effects. Of course, the gluino exchange diagrams proportional to a? tlong ma y 
be particularly relevant. There are several sources of flavour mixing in gaugino couplings that we 
now turn to discuss. 

Within the general framework of supergravity [15], a theory is defined by the gauge and mat- 
ter superfields, by their couplings to supergravity encoded in the Kahler potential, and by the 
generalised gauge and Yukawa couplings. The low-energy theory is then fixed by the values of 
the auxiliary fields which provide the spontaneous breaking of supersymmetry. Therefore, the 
parameters in the low-energy theory with broken supersymmetry are to be determined within the 
framework of a larger theory explaining all interactions including gravity. Superstrings [16] are the 
best prototypes, despite being incomplete. Nevertheless, recent studies [17 — 19, 12] of the field 
theoretical limit of superstrings have given hints on the general form of the relevant parts of the 
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Kahler potential as well as on their interpretation in terms of stringy symmetries. A conspicuous 
result of superstring model calculations is that the three families of quark superfields may couple 
to supergravity according to quite different terms in the Kahler potential. Therefore, the minimal 
hypothesis in phenomenological analyses, alleging that supergravity couplings are flavour indepen- 
dent is not quite confirmed by these explorations within superstring theories. The dilaton superfield 
in these theories does have universal supergravity couplings to quark superfields [8,20]. But the 
moduli fields associated to the compact 6-dimensional manifold, needed to reduce the space-time 
dimension to four, have model dependent couplings to the quark superfields [17, 18, 12]. Thus, the 
Kahler potential can be different for each flavour. Moreover, gauge singlet fields have been conjec- 
tured to explain the so-called textures in the Yukawa couplings [21,22]. They are assumed to get 
relatively large vacuum expectation values and could introduce non-universality corrections to the 
supergravity couplings of quark superfields. 

At the level of the effective (renormalisable) theory, below the Planck scale, the supersymmetry 
breaking effects reduce to gaugino masses and the soft interactions in the scalar potential [23]. There 
appear hermitian (mass) 2 terms involving the scalar fields of each given SU(3) x SU(2) x U(l) 
quantum numbers and their complex conjugates. The corresponding sub-matrix of the scalar 
(mass) 2 matrix depend on the Kahler potential and on the supersymmetry breaking auxiliary fields 
[24, 15]. In particular, the universality or flavour independence hypothesis assumes equal masses for 
all squarks at the unification. At lower energies, this conjectured universality is broken by radiative 
corrections due to Yukawa interactions and acquire some calculable flavour dependence. Moreover, 
analytic (mass) 2 terms involving scalar fields of opposite SU(3) c xU(l) em quantum numbers are 
generated after the electro-weak gauge symmetry breaking by the Higgs scalars. They are basically 
proportional to the Yukawa couplings in the superpotential. Again, if universality is assumed for 
the proportionality factors, referred to as ^-parameters, at the unification scale, their equality is 
spoilt at lower energies by the calculable radiative corrections. 

In most discussions of supersymmetric FCNC processes, universality of soft terms is assumed. 
Then, the most striking effects of radiative corrections are of two kinds: a) gauge corrections which 
are universal and, in general, tend to attenuate loop effects by an overall rise in the squark masses 
if gauginos are relatively heavy, and b) Yukawa corrections dominated by the top coupling, A t , tend 
to align the down squark mass eigenstates to the up quarks [3 — 5,7]. This reverses the pattern 
of gaugino couplings in comparison with the gauge boson ones. Chargino couplings to down 
squark and up quarks are approximately family diagonal while gluino and neutralino couplings 
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become proportional to the Kobayashi-Maskawa matrix. However, the expected physical effects 
are generally quite tiny with the present overall bounds on supersymmetric particles 1 . The b — ► 
57 transition is the remarkable exception which gives interesting information on supersymmetric 
theories [7,8,10]. 

Motivated by superstrings, as well as symmetries proposed to explain the structure of Yukawa 
couplings, new analyses [11 — 14] have been performed on FCNC transitions produced by non- 
universality in supergravity couplings. If the flavour dependence is large already at the supergravity 
or Planck scale, the quark mass eigenstates and the squark mass eigenstates are not superpartners 
under the action of supersymmetry. There is a unitary transformation in the family space reflecting 
the difference of the basis for physical quark and squark states. The gluino and neutralino couplings 
will be proportional to this unitary matrix. The natural supergravity basis for quarks and squarks 
is established by diagonalisation of the Kahler metrics (with its field-dependence classically fixed). 
In this supergravity basis, the squark (mass) matrix is dominantly diagonal but the eigenvalues 
may be different. In the same basis, the Yukawa couplings are expected to be non-diagonal. In 
order to define the quark mass eigenstates, one has to perform four unitary transformations to 
diagonalise the up and down Yukawa couplings in the family space at the supergravity scale. 
These unitary matrices tend to produce FCNC effects in neutral gaugino exchange already at tree 
level. But the difference of many orders of magnitude between the unification scale and the low 
energy scale suitable for FCNC calculations, suggests that the quantum running of the quark and 
squark masses and mixing matrices could affect the tree-level pattern. 

In this paper, we present a detailed study of the renormalisation group equations for the 
Yukawa and soft-terms in the quark-squark sector, written as general matrices in family space [4, 5]. 
The phenomenological hierarchy in the Yukawa coupling eigenvalues is exploited using the obvious 
device of keeping the up quark Yukawa couplings diagonal through scale dependent rotations. 
The RGE for Yukawa couplings and soft terms are then easily solved in the approximation in 
which all Yukawa couplings but At are neglected. It is straightforward to compute the corrections 
proportional to the other couplings (A^,A^, etc ) to the lowest order. Although they may be 
relevant to some applications, these corrections are not included in this paper. 

In Section 2, the RGE for the Kobayashi-Maskawa matrix is derived, then solved in the case 
in which only the third family Yukawa couplings, A t and Xb, are retained in the RGE. These results 
already exist in the literature [4,25], but here, the derivation is more transparent. In particular, 



1 For a recent discussion, see e.g., Refs.[10, 11]. 
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the role of the Yukawa coupling hierarchy is exhibited so that the accuracy of the approximation 
can be readily estimated. 

In Section 3, the RGE for the squark soft analytic terms are analysed. We argue that the 
current way [15] of expressing their proportionality to the Yukawa couplings should be changed 
to better implement the general boundary conditions at the unification or Planck scales. Our 
parametrisation is more suitable because Yukawa couplings relate left-handed and right-handed 
quarks whose supergravity properties may be very different. The RGE for the newly defined 
matrices of parameters are displayed and analytically solved, in the limit of A t dominance, in full 
complex matrix form. 

In Section 4, the RGE for the hermitian sector of the (mass) 2 matrices are analytically solved 
within the same approximation. 

The calculation of the gluino loop contribution to the electric dipole moment [26] of the neutron 
is revisited in Section 5, in view of possible effects of non-universality in the soft terms. We notice 
that contributions proportional to m& or m t which are usually neglected, could overcome current 
contributions proportional to ma or m u . In order to estimate their relative strength, the upper 
limits [8, 10, 11] on the flavour mixing components of the scalar (mass) 2 matrices are used. 

In Section 6, we explore the effect of the RGE evolution on FCNC effects in the presence of 
flavour dependent supergravity couplings. In the absence of a realistic model, we define an ad hoc 
one, where the mixing matrices relating the basis that diagonalise the Kahlerian metrics and the 
Yukawa couplings, respectively, are maximal. With simple assumptions on the Kahler potential 
and the superpotential, we find that in some extreme cases the existent upper bounds on FCNC 
generating parameters can be saturated. 

2. QUARK MIXING MATRICES 

Let us begin with the RGE for Yukawas and KM matrices. We define the unitary matrices 
Ur, Ul, Vr and Vl by the diagonalisation of the up and down Yukawa couplings: 



where all the parameters have their dependence on t = In (A/A ) / (47r) 2 , defined by the well-known 
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(2.1) 



RGE 2 



dt 
dX D 

dt 



3 (a^ + Tr + X D X D - 2C u a g 2 

d [3 (x d X d + Tr A^A D ) + \\j\u + Tr A^A L - 2C d a g 2 a 



A, 



16 9 9 13 9 
y 93 + 3^2 + 

16 o o 2 o 



(2.2) 



It is convenient to write \u and Ad in terms of the hermitian matrices 



Yu — XjjXu , Yd — X d Xd 



(2.3) 



diagonalised by the left-handed unitary transformations Ul and Vl only, and XuX\j, XdX d diago- 
nalised by Ur and Vr. We then work with (2.3) only, the right-handed counterpart being analogous. 
The RGE can then be cast in the form: 



dYu 
dt 

&u 
dY D 
dt 



{3Y U +Y D ,Y U } + 2A U Y U 
3TtYu- 2C u a g 2 a 
{3Y D + Yu,U D } + 2A D Y D 
3 Tr Y D + Tr Y L - 2C°g 2 a 



(2.4) 



Now diagonalise Yjy and Yd : 

% = u}YuU L 



Y d = VIY D V L = 

and define the KM matrix (an analogous one may be defined for the right-handed states) 

K L = VlU L 

Then the RG evolution for the eigenvalues in Yjy follows from (2.2): 



A? 



X 2 



X\ 



A? 



(2.5) 



\2 



(2.6) 



dY v 
dt 



6Y V + 2A V + 2 



(k{y d 



K, 



(d) 



Yj 



u 



(2.7) 



2 See, e.g., Ref.[4]. 
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where we denote the separation of any matrix into diagonal and non-diagonal parts, M = + 
M( n d)- Consider, e.g., U^^Ul, which is an element of SU(3)/iT where H is the SU(3) Cartan 
algebra. The off-diagonal counterpart of (2.7) coming from (2.2) reads: 



(k{y d k l ) % 

\ / (na) 



(2.8a) 



Notice that K\YdKl is just Yd in the basis where Yjj and SU(2) wea k are diagonal (which will play 
a major role in the next sections) at each t. Analogously, 



V, 



dV L 
dt 



,Y D 



K[YuK L ) Y D 

(na) 



Finally, 



(2.86) 



(2.9) 



dt L dt 

At this point we can take advantage of the strong hierarchy of the (mass) 2 eigenvalues for up- 
quarks on the one hand (m 2 » m 2 » to 2 ) , and for down-quarks on the other hand (m 2 » m 2 S> m^) 
Taking the components of (2.8) 



U 



t dU L 



V, 



dt 



dt 



ab 



m 2 + m 2 
m 2 — m 2 



a 6 



(a, b = u,c,t; a ^ b) 



ab 



m\ + m 2 

m 2 — ml I \ ~ "/a6 



K L Y D Kl 



(2.10) 



(a, b = d,s,b; a ^ b) 
one obtains the evolution of l7 l and Vl from the (mass) 2 hierarchy: 



dt 

L dt 



a 6 



\ / ab 

U (k l YuKI) 

\ / ab 



(2.11) 



(m a < m 6 ) 
(a = 6) 
(m a > m 6 ) 



Next order approximations in (m 2 /m 2 ) and (m 2 /m 2 ) are easily obtained (but more difficult 
to integrate), though useless in physical applications. From (2.9) and (2.11), we derive the RGE for 
the elements of the evolving KM matrix. We fully use the hierarchy among masses of equal-charge 



quarks and get the following results 3 4 : 

d 



dt lnK bt = (+(\ 2 t +\l) 
±\nK bc = ( + X 2 b \K bu \ 2 
^-lnK st = ( + \ 2 t \K dt \ 2 

dt s * 1 1 (2.12) 

±lnK bu =( + \ 2 b \K bc \ 2 

j t lnK dt =( + \ 2 t \K st \ 2 

C = -(X 2 + X 2 )\K bt \ 2 

Unitarity is easily checked. The RGE for remaining elements of the KM matrix are more 
involved, so we take advantage of the strong hierarchy in the non-diagonal elements to write them 
as follows: 

—In K du ~ A^ \K SU \ 2 + X 2 b \K bu \ 2 

+ X 2 c \K dc \ 2 + X 2 \K dt \ 2 
^\nK sc ~X 2 s \K su \ 2 + X 2 \K bc \ 2 

+ Xl\K dc \ 2 + X 2 \K st \ 2 (2i3) 
-\nK dc ^-{X 2 s +X 2 c )\K sc \ 2 



+ {X 2 + X 2 )\K dt \ 



2 



^\nK su ~-{X 2 s + X 2 c )\K sc \ 2 
+ {X 2 b + X 2 t )\K bu \ 2 

These results have been obtained by a different method in Ref.[25] 5 . It is notice- worthy that 
the present approach clearly shows that the relevant approximation is given by (2.11). The RG 
evolution of the KM matrix elements in (2.13) is negligible compared to those in (2.12). The latter 
are easily integrated if one neglects the small matrix elements in the r.h.s. . In practice, it is enough 
to put \K bt \ 2 ~ 1 to get 

K bu (A) K bc (A) K st (A) 
K bu (Ao) K bc (A ) K st (Ao) 

e It+I » (2.14) 



K dt (A) _ Jt+Ib 



K dt (A ) 

It = J dt X 2 {t) h = J dt X 2 (t) 



3 Notice that in spite of their real aspect, Eqs.(2.12) and (2.13) are RGE for complex quantities. 

4 For simplicity the index L is omitted. 



5 And, previously, in Ref.[4] for X 2 » A^. 



while all other KM matrix elements have a negligible ^dependence. Assume Af S> \%, then, from 
the RGE for Af , one gets 



(6 3 = 3, b 2 = -1, h = -11) 

where Aj? rit (i) is the A^ value corresponding to a Landau pole at to. Hence, the RG evolution of the 
KM matrix can be dismissed in most problems, or taken into account by simple approximations 
as in (2.14). 

It goes without saying that everything can be repeated ipsis litteris for the unitary transfor- 
mations Ur, Vr, Kr on the right-handed quarks. 



3. ANALYTIC SOFT TERMS 

In this Section, we establish the RGE for the so-called soft analytic scalar interactions. Those 
relevant for the phenomenology of flavour mixing are: 

H^rjuQ + H 2 D C VD Q + h.c. , (3.1) 

where isospin indices are omitted, r\u and r\u are matrices in the family indices. The matrix el- 
ements of i]u and t]d are proportional to the supersymmetry breaking scale, hence to m 3 / 2 . As 
discussed in the introduction (see also (3.5) below), it is simple to assume them to be also pro- 
portional to the corresponding matrix elements of the Yukawa couplings, A^ and A^, respectively. 
The current assumption of "universality", i.e., flavour independence of the soft terms is usually 
stated by setting [15] i]u = A\u, r]o = AAd where A is (m 3 / 2 ) ■ Since flavour independence 
is broken in the RGE by the effect of Yukawa interactions, this is usually taken into account by 
replacing A by a diagonal matrix and neglecting the scale dependence of the KM matrix. In more 
specific calculations (e.g., in that of the electric dipole moment of the neutron in models where 
A is assumed to be universal) the matrix character of the soft terms have to be more carefully 
manipulated [27]. Nevertheless, the soft couplings i]u and r\u are not expected to be in general ex- 
actly flavour independent in supergravity theories. In such circumstances, one should choose some 
parametrisation of r\u and r\n which simply embodies the assumed dependence on \y and Ad and 
still preserves the generality of supergravity breaking terms. Obviously, it is always possible to 
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define a matrix A as A 1 rj (or 77 A 1 ), but since these soft terms and the Yukawa couplings link 
quarks with different ^-numbers (L and R states, respectively) let us scrutinise the RGE for r) u , 

^= V U [5X U ^X U + A Dt A D + 3 Tr (X^X U ) - 2C^l] ^ 
+ 2X U [2X U ^ U + 3 Tr (X u ^ u ) + A Dt r/ D - 2C^g 2 a M a ] 

and its analogous for t]d, replacing U ^ D (and adding the leptonic couplings Tr A L ^A L ). Actually, 
the choice rj u = X u A u and r\ D = X D A D is quite a satisfactory definition leading, e.g., to the 
following equation for the matrix A u 

(\A U 

_ = 5A U Y V + YuA u + Y D A D + 6 Tr {A U Y V ) - ^g 2 a M a (3.3) 

This choice has been used in the literature 6 in connection with FCNC effects. 

However, in the more general context of flavour dependent supergravity couplings, the choice 
in (3.3) is not the best in view of the structure of the boundary conditions at the unification scale 
A , not far from the Planck mass. Indeed, consider an effective supergravity theory [15] defined 
by some Kahler potential K (X, Hi, H 2 ,Ui,Qi, Di„.) (i = 1,2,3 is the family index) and a 
superpotential: 

W = H 1 U C X U (X)Q + H 2 D C X D (X)Q + ... (3.4) 

where we have defined the chiral field X as the goldstino (i.e., supersymmetry breaking) direction, 
such that the auxiliary field Fx ~ (m 3 / 2 ) (in units Mpi anc k = 1). Then, the trilinear coupling 
matrix rj u will take the form [15, 24] : 

^ = {(h + \m) ^ {x) + r ™W x) + r « A & w + r ^ (*)} f x 



l xh~ k ~dX 



(3.5) 



where Khh*, K(.*i and Kg*i denote the effective Kahler metrics for Hi, U and Q, respectively. 
This formula clearly shows how flavour mixing can be introduced in the 17's, through either the 
"Kahler connexions" T,T or a non-trivial dependence of the Yukawa couplings on the goldstino 
scalar partner, X. 



See, e.g., Refs. [4, 8]. 
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Therefore, we favour a different decomposition of the soft terms, as follows, 

V u = B U R \ V + XuBY 



(3.6) 



whose advantages become clear comparing with (3.5). This decomposition is defined up to terms 
commuting with \jj (Xd, resp.) which can be arbitrarily transferred from Brs to B^s, or vice- 
versa. In our expressions there below, we shall equally share the family independent parts between 
Br and Bl- The RGE for the Br's and B^s immediately follow from (3.2). However, it is more 
convenient to treat these RGE in such a way as to easily utilise the hierarchy of Yukawa couplings, 
in particular, the large value of the top mass. The obvious idea is to keep A;y(A) diagonal at any 
scale A by the transformations Ur, Ul defined in (2.1). This amounts to choosing a rotating basis 
because of (2.8a) and its analogue for Ur. Let us denote the diagonal Yukawa coupling matrix by 
Afy(A) = U R (A) \u (A)Ul(A) , and any matrix in flavour space in the corresponding basis by 



M(A) = £/t(A)M(A)£/(A) 



(3.7) 



where U = Ur or Ul, according as M acts on R or L states, respectively. Hence, one derives the 
RGE for M(A) in the rotating basis as 



dM 
~dt~ 



C*(A)^E/(A) + 



(3.8) 



where the last term takes into account the rotating basis. From the discussion of the previous 
Section, one knows that the effect of this last term is quite small in the L-sector, where the KM 
matrix in (2.11) is known. On the contrary, the corresponding KM matrix for the i?-sector is 
unknown and so are the angular velocities defined by the analogue of (2.11). 

Let us now write down the RGE for the B's in (3.6). Unless otherwise stated, all matrices are 
written in the rotating basis (diagonal A^). (To simplify the notation we omit the " " symbol on 
these matrices). From (3.6), (2.2) and (3.2) one then derives the following RGE equations in the 
(real) diagonal \u basis: 



dB 



u 



dt 
dBY 
dt 



Z = 2B u R \l + A\lB u R + 6 Tr (A 2 ^) - 2C^g 2 a M a t + 



t dU R 



dt 



bB u L \l + \lB u L + 6 Tr (\lB u L ) - 2C u a g 2 a M a t + \b%, A^aJ + 2\^ D B°\ D (3.9) 



+ 2\ d \dBl + 



dU L 
■ dt 
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2B%\i D \ D + 4\ i D \ D B% + 



D 



dBl 
dt 



dB ° - 5BE\1;X d + \\j\oBE + 



6 Tr ( X^XdB^ 



2C d a g 2 a M a 



1 + 



°Ri U R 



dt 



dt 



6 Tr {X^ D X D BE ) 2C d a glM a ] 1 + [Bf , A^] + 2Agi?gA l/ 



+ 2A 2 ^ + 



dZ7 L 



where M Q are the gaugino masses and the terms proportional to the identity matrix 1 have been 
equitably distributed between Bl and Br. The last term in each equation corresponds to the action 
of an element of the SU(3) algebra whose expression is given by (2.11). Notice that U ] L (dU L /dt) 
is known in terms of the KM matrix but U R (dC/^/di) is unknown. 

Let us now illustrate some of the feature of (3.9) by a simple, yet possibly realistic, case where 
the RGE can be solved analytically. This is so if all Yukawa couplings can be neglected but the 
top one, A t . This approximation which kills many interesting terms in the RGE above, allows one 
to compare with previous work 7 in the literature where family mixing is neglected. Indeed, the 
solutions of (3.9) in this approximation are as follows: 

{B U L + B U R ).. (t) = {B U L + B U R ).. (0) - 4 £ CZ&tMo 

a 

-§(l + *a)[(B£ + BK)33 (°)+^o] 

(sKL(f) = e-°8* (BKUO) 



(B U L ) ij (t)=e-^< (0) 



"I 



"32 



"f 3 



"fl 



2 
4 



"23 



"32 



"f 3 



"31 



a 



R 
12 



a. 



21 



a; 



12 



a; 



21 



(Bj?) y (t) = . . (0) - My £ C^Mo 



(i,j = 1,2,3) 

(Sf ) . . (t) = (B?).. (0) - 2 J] C&frMo " ^Us [(BK + ^) 33 (0) + CM ] 



(5f) 12 (*) = (Bf) 12 (0) 



B 



' i3 



(t) 



-It 



B 



' iZ 



(0) 



1,2 



(3.10) 



(BE) 3l (t) = e I *[(BE) 3i (0)-2l(BY) 3i (0)] , 



7 Analytic solutions without family mixing in the A t -dominance approximation can be found in 
[4] and [28] with universal supergravity couplings. 
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where I t and e Tt are defined in (2.15), M is the gaugino mass at the unification scale (t = 0), 
and p and £ are defined by: 

p = = 1 _ e" 12 ^ 

AcritW (3-11) 

£ = -12A 2 rit t-l . 

Notice that e 7 * ~ 1, so that the sizeable effect in the B£ , B^ matrix elements are the gauge 
contributions in the diagonal and the top trilinear term in (5^ ) . 

The consistently complex character of all the relations in this section is notice-worthy. 

4. SCALAR MASS MATRICES 

The RGE for the scalar (mass) 2 soft terms are available in the literature in full matrix form 
[4,8]. However the universality hypothesis has been used in the solution [28]. The one- loop RGE 
for the relevant masses are, 

— -9. = |a^A[/ + \ ] D \ D ,m 2 Q^ + 2 (m 2 Hl \\j\u + m 2 H ^ D X D + X^mljXu + X ] D m 2 D X D 

WuVu + vivo) ~ 8C%g 2 a M 2 
^ = 2 {XuX^ml} + 4 (m 2 Hl XuX\j + Xjjm 2 Q X\j + VuV l) - 8C^ 2 M 2 + ... (4.1) 

= 6Tr \ rn 2 Hl X\jXu + m 2 Q X^Xu + m^A*, + r^rju) - 8C^g 2 a M 2 + ... 

where Tr is the trace on flavour indices. The RGE for m 2 D and m 2 H2 are obtained from those of 
mfj and rn 2 H , respectively by the exchanges U < — ► D, H\ < — ► H2 and inclusion of the leptonic 
terms. Because of non-universality, there is the additional term YiS in the RGE, where, 

S(t) = I^Tr (Ym 2 ) (4.2) 

The last term in each RGE takes into account the change of the basis at each energy if, for instance, 
Xjj is kept diagonal. Let us then display the solutions of the above RGE under the assumption 
that all Yukawa couplings but A 2 are neglected, consistently with (3.10). In the basis in which Xu 
is diagonal, we obtain for those hermitian matrices: 

2 

,2\ n\ _ ( m 2\ (n \ j_ 9i+ Q(t\\ _ P k _ ^,Qi 



K) 33 (*) = K) 33 (°) + T s (°) ~i K - * * M o 

|^5(0)-f 



K) 33 (t) = (m 2 / ) 33 (0) - \g\t 5(0) - P -K - 7 u t M 2 
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ml(t) = m 2 D (0) + (M t S(0) - 7 D t M 2 ) 1 

(ml)., (t) = (ml).. (0) + &t S(0)5y + (e" 2 ^ - l) (^) 3j (0) - M 2 b %] 

(ml)., (t) = (ml).. (0) - -g\t 5(0)% + (e" 4 ^ - l) (B u R ) iZ (B U R *) j3 (0) - 7 u t M 2 8 %3 



K) l3 (*) = e " 7t K) l3 (o) - K% (o) 



p 

K) i3 (t) = (ml) i3 (0) - P -Al (B%) a (0) (4.3) 
where i,j = 1,2. We have defined the following quantities: 
A = (B U R +B U L )^ (0)+£M o 

K = (m 2 Q + ml + m 2 H2 l) 33 (0) + (1 - p)\A \ 2 - £ 2 M 2 



-8^o 2 (e + i)tM ( 



7 Q = y<?l (1 + 3</|t) + 6 52 2 (l - <&) + \dl (l - Us?*) 
7 17 = f </ 3 2 (l + 3s§f) + f g\ (l-llrft) 

7 D = f 5 3 2 (l + 3^) + ^ 2 (l-ll^) (4.4) 

As already discussed in the literature [3,4,8], starting from universal soft terms, the flavour 
dependence due to the RG running only appears in the diagonal matrix elements. In the general 
case instead, the flavour dependence of the (mass) 2 terms gets contributions from the BY, B R 
matrix elements. Also in this Section, the complex nature of the various parameters is explicitly 
respected. 

Let us estimate the impact of the running of the soft terms from the unification scale down 
to the SU(2) x U(l) symmetry breaking scale. Assuming h 2 ~ 1 at A ~ 0(Mz) , one gets the 
numerical values: 

t = -.2 p = .8 £ = 2 
e~ h = .87 I t = .13 

2j2c-g 2 J = 2j2C d a g 2 a t = -1.9 
7^ = 6.7 7^ = 7^ = 6.3 (4.5) 
The substitution of these numbers in (3.11), (4.3) and (4.4) gives the following relations, 
(BY + B U R ).. (t) = (BY + B U R ).. (0) + 4M - .4 (1 + 5 t3 ) A 
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(Bg + B£) u (t) = (BE + Bg).. (0) + 4M - .135 i3 A (i = 1, 2, 3) 
K), 3 W-5K) l3 (0) (i = l,2) 
W) 33 (*) = ( m Q) 33 (°) - - 4 ™o " -025^ + 5.7M 2 
(mj) . . (t) = (m 2 Q ).. (0) + 6.7M %- - .25 (#£% (i#) 3 . (0) 

K) l3 (*) = - 87 K?) i3 (°) - - 13 ^o (BY% (0) = 1, 2) 

W) 33 (*) = W) 33 (°) - - 8rn ° - mA ° + 43M o ( 4 - 6 ) 

K) y (t) = (m 2 ,).. (0) + 6.3M % - .4 (Bg) . 8 

(m 2 / ) i3 (t) = .87(m 2 Q ) i3 (0)-.13A* (^(O) 



where, 



m o = \ ( m Q + ™u + m 2 H1 ) 33 (0) 
A = (B U L +B U R ) ZZ (0) + 2M 

and Mo is the mass of the gauginos at the unification scale, assumed to be all equal. 

Since FCNC effects from the non-universality of mfj, m 2 D , m 2 Q are expected to be relatively 
small, it is convenient to expand the contributions to the possible FCNC processes in powers of 
the following matrices: 











(*^r) ij 
















(4.7) 



where 

m 2 Q = ^Tr m 2 Q (0) + "f Q \t\M 2 + 9 X S(0) 

ml = ^Tr ml(0) +j u \t\M 2 - ^g 2 t S(0) 

(i.e., only the gauge contribution in RGE are included in the definition of these "average" masses). 
By definition, these matrices are defined in the diagonal Xu basis. Instead, for the corresponding 
5r and 5£, it is more convenient to go to the diagonal Ad basis through KM rotations, Kl, Kr. 
Hence, we define the following matrices: 

(.K R m 2 D K R ) 
V R )ij~ m 2 D 6 13 

(K L m 2 Q Kl) 

\ / ij mQ 
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where the "~" symbol indicates the diagonal Ad basis for the matrices. Unfortunately, 5 R de- 
pends on the unknown K R matrix. The matrices defined in (4.6) are those usually compared to 
experimental data on the various FCNC processes. Thus, there are several upper bounds in the 
literature obtained by considering a variety of physical transitions [8, 10, 11]. However these upper 
bounds correspond to different assumptions and approximations, and also depend on the overall 
squark and gaugino masses. 



5. AN EXAMPLE: THE NEUTRON E.D.M. 



The neutron electric dipole moment is a precious parameter in fundamental physics as it is 
strongly suppressed in the Standard Model, on the contrary, it is sensitive to new sources of CP 
violations. Among the possible physical phases in supersymmetric theories, some are more relevant 
for the discussion of the neutron e.d.m., namely, the relative phases between the matrix elements in 
the soft analytic scalar couplings and the gaugino masses [26, 27]. We choose the usual conventions 
and redefine the fields so that gaugino masses are real. In this Section, we choose to work in the 
basis where the down-quark mass matrix is diagonal and real. With respect to the definition in 
the previous Sections, we have to perform KM rotations to define, e.g., 

m 2 Q = K L m 2 Q K{ 
m 2 D = K R m 2 D K R 

VD = K RVD Kl (5.1) 

Since the calculation of the gluino-squark contribution to the e.d.m. D n of the neutron has 
been already presented many times in the literature, we just quote the simple final expression for 
the e.d.m. of the d-quark: 

D d = ~e^MiV 2 lm T dd 

f k 2 dk 2 ( 1 1 \ 

Ud -J [k 2 + M 2 }{ki+rn 2 DR VD ki + rh 2 D J dd ^ 

where M3 is the gluino mass, only the lowest order in the analytic trilinear coupling is retained, 

the phase in the "/i — term" is neglected for simplicity, and the matrices m 2 R and m\ are defined 

as: 

2 M 2 



2 ~ 2 

1 Dr =™ 



2 d+(XdV 2 ) - sin 2 6 W cos 2 



m 2 DL = m 2 Q + (\ D V 2 ) 2 - M% Q - cos 2 /3 (5.3) 
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tan (3 = , Vi = (Hi) 

V2 



(All matrices are denned in the diagonal Ad basis). 

In order to extract the dominant contributions in the presence of flavour dependent super- 
gravity couplings, let us consider a simple approximation. We consider a Taylor expansion of the 
propagators (k 2 — m 2 DR ) 1 and (k 2 — rh 2 DL ) 1 around some mean value m 2 , an average over the 
squark and gluino masses, and then we integrate on k 2 . By keeping only the dominant term of 
each kind, using the notation introduced in (4.6), one ends up with the following expression: 



(5.4) 



■— l - Im 
30 



where i = b,s; j = b, s, d. Because of the Yukawa coupling hierarchy, all terms in (5.4) can be 
equally relevant. Indeed, let us consider the upper bounds obtained in Refs.[8, 10, 11] for the 
squark mass differences, assuming fh ~ 0(1 TeV). They are as follows: 



(if) 



sd 



< .1 



< .3 



°R 



bd 



ds 



< .1 



ds 



< .3 



sd 



(*«L(^) 



bd 



db 

< 4 x 10 -5 

< 5 x 10" 3 



(5.5) 



Combining these values with the quark mass ratios, we get the following estimates for the relative 
magnitude of the various terms appearing in (5.4) 8 : 

12 m s 



20 m d 
12 mi, 

20 m d 



12 m s 

30 ma 
12 m b 

30 m d 

Though the coefficients of [B£ 



(*f(fl)) 



ds 



db 



0(1) (0(-07)) 
0(70) (0(20)) 



~ (3 x 10" 4 ) 
are allowed to be rather large, upper bounds 



BR 



are also available for these matrix elements. A rough upper bound is obtained by assuming 
8 The estimates in brackets are obtained from the bounds on the products S^SR and the 



assumption 8® ~ 5 



D 
R- 
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[B^ R ^J ~ (in) , where m is the average mass in (5.4) and then imposing bounds on (n)bd, (v)db 
obtained from studies of supersymmetric effects in K and B physics [8, 10, 11]. However, once trans- 
lated in terms of the parameters in (5.4), these bounds just require (Bl{r) ) ~ ^ (j3£ itj ) • 
Thus, terms proportional to nib can even dominate the first, more familiar, term. 

We have dealt with Dj, the e.d.m. of the (f-quark as some comparison of the different contri- 
butions can be performed. However, the importance of terms proportional to heavy quark Yukawa 
couplings is even more obvious if one calculates D u , the e.d.m. of the u-quark. In this case we 
work directly in the diagonal \jj basis used in the previous Sections. The analogue of (5.2) reads 

1 

D u = +—e^-M*V 1 Im T uu 
f k 2 dk 2 ( 1 



[k 2 + M 2 f \k 2 + m 2 UR k 2 + m 2 UL 
VI 
3 



IM 2 

m u R = m u + {^uVif + —Z- sin 2 6 W cos 2(3 



m 2 u L = m 2 Q + (XuV,) 2 + M 2 Q - cos 2 (3 (5.7) 

Then the same approximations as in (5.4) (but, in this case, they could be poorer) yield, 

* Im ^ ± Im (BK + Bl) uu - ^ Im [(^ + 

(B U L ) l3 + (<5K) U , (Sl) iu } (5.8) 



mi 

Im 

30 



Estimates like those in (5.6) are only possible for (uc) or (cu) matrix elements. Using results 
from [8, 10, 11] we obtain the very poor constraints, 

12 m r 



20m, 



~0(14) (0(6)) 

,Ulli u \ 

12m 



it. C 



K) uc (S U l) c -0(.2) (5.9) 



30m ' " "'• ' 



The largest contributions are certainly those proportional to m t , the top mass. Let us scrutinise 
the RGE of the relevant soft parameters in the approximation of (3.10). First, we notice that 
(.B^) t an d C^l)* are only slightly renormalised even for Aj near the critical ("fixed-point") value 
(while the transposed matrix elements get more reduced). Instead the (B^ + B^) (33) diagonal 
element can be strongly reduced in the large At limit so that non- diagonal and (BjQ can 

dominate. This is even enhanced by the large contributions to (5^) , (<5£) diagonal terms for large 
A 2 - 

Therefore, the contributions to Dd and, especially, D u from terms proportional to m& and 
m t , respectively, could be larger than those obtained by neglecting the squark mass splitting in 
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the propagator, which are proportional to light quark Yukawas. However, the whole discussion 
remains merely academic in the absence of a specific model for supergravity parameters, including 
the phases in the soft terms. The minimalist assumption is to require universal supergravity 
boundary conditions and the KM phase as the only source of CP violation. This has been shown 
to yield extremely small contributions to the quark e.d.m. [27]. Another rough estimate [26] consists 
in setting (Br + Bl) ~ M3 ~ m to obtain for the phase <fi of (Br + Bl) , sin (\> < (700 MeV/ra) 2 , 
from the experimental limit on Dn- These two extreme cases clearly show that the effects could 
be important compared to the Standard Model prediction where D N is very strongly suppressed. 

6. MAXIMUM FLAVOUR MIXING IN SUPERGRAVITY 



In this Section we aim at studying a model with maximal mixing at the supergravity scale, 
in order to get some quantitative insight about the effects at low energies. A natural model 
for flavour mixing would constrain both the Kahler potential and the superpotential (Yukawa 
couplings) in terms of symmetries. We do not know a realistic model of this sort and so we forget 
the naturalness requirement. Thus, we discuss a model designed to stress FCNC effects where 
the flavour dependence is introduced by hand in the (field dependent) supergravity metrics. For 
simplicity, the Yukawa couplings are assumed to be left-right symmetric. The U and V matrices 
at the unification scale Agut are defined from the Yukawa couplings at Agut by (2.1) and the 
KM matrix is given by the relation 



V = U K\ 



(6.1) 



Let us consider the real matrix 



y/6 



( Vs 1 V2\ 

-y/3 1 V2 (6.2) 
V -2 y/2/ 

(we neglect all phases in this Section). From (6.1) one gets the following model for \u (Agut) : 

/ 1 + e u + e c l-e u + e c 1 - 2e c \ 



h 
3 

3A, 



2 A 



1 E u H - £ c 

V l-2e c 

Sr. = 



l + e u + e c 
l-2e r 



1 - 2e r 



1 + AeJ 



2X t 



(6.3) 



19 



This matrix is (almost) "democratic", in the sense that all matrix elements are almost equal. In 
this sense the flavour mixing is maximal. Then for the sake of our qualitative discussion , we define 
V as in (6.1) with K phenomenologically given. 

The supergravity basis is defined by the diagonalization of the Kahler matrix. In this basis, 
an honest assumption is to set all soft terms diagonal at Aqut (however, this is not necessarily so, 
as discussed in Ref.[6]). Thus, let us define 



(BY) 



13 



A u 



+ b 



u 



Li I "1] 

( m u)ab= {™U + VUa)S a b 
( m D)ab = ("& + A*L) S "b 

(i,j,a,b= 1,2,3) 



(6.4) 



By performing the U rotation, the soft terms are transformed into the more convenient diagonal 

u 

R 



Xu basis. In this basis B R and B^ take the form: 



/ V3Ab u 



B U R 



Ab R 
A% 



2 3 



u 

Ra 



V2 



R 



A% 



\ V2A'b u R 



b U R1 +b u R2 
b u -h u 

°R1 °R2 



2^3 



A% 



V3A6g 



Ab R 



V2A>bU \ 



Ab R 



/ 



(6.5) 



and the analogous expression for B% with R — ► L. The corresponding relations for B R , B® are 
analogous with U — > D in the basis where is diagonal and thus they are transformed by a 
KM rotation into the diagonal Xjj basis. 

Similar expressions also result for the m 2 parameters, which can be written as: 



(>-»? + i£ 



( V3A/ij 



\V2A'- 2 



(I = Q,U,D) 



Atf 



V2A'tf \ 
Atf 
/ 



(6.6) 
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In order to proceed and fix some parameters in view of a quantitative analysis, we consider a 
simple Ansatz for the quark quadratic terms in the Kahler potential 

where $ A = Q l , U a ,D a . The moduli fields introduced in the effective supergravity by the compact- 
ification of six dimensions in string theory, are represented here by the overall modulus, denoted 
T. The exponents tia are the modular weights that specify the behaviour of each superfield $a 
under modular transformations of the compact manifold. We shall assume them to be negative 
integer numbers [17, 18]. Then, from the general supergravity expressions: 

m AA* = (F M F N *) (Zmn*Zaa* — Rmn*aa*) (6.8) 



Vabc — F D m Xabc 

Dm = Z- 1 (d M + \d M K^j Z 
Zab* = d A d B *K (6.9) 

one gets the soft scalar couplings and masses in terms of the supersymmetry breaking auxiliary 
fields F M . Here we consider the simplest case in which these auxiliary fields are just F s and F T , 
corresponding to the dilaton sector and to the overall modulus sector in the effective supergravity 
theory obtained from the orbifold compactification of superstrings. Accordingly, we introduce the 
Goldstino angle 6 [12], such that tan# = F s /F T , and neglect possible phases. Then the gaugino 
masses at Aqut are all equal and given by Mq = | F s | 2 , while for the other parameters one obtains 
[12] 

M 2 

m\ = (1 + (1 + n A ) cot 2 6) 
{B u R ) ab = -M,(\ + + + n u a ) cot 9)5 ab (6.10) 

with analogous expressions for B^, Bf. 

These expressions can be substituted in (4.6) to allow a numerical investigation. The integers 
rii (modular weights) are free parameters. We choose to fix their values n\ = —1, 112 = —2, = —3 
for the first, second and third families. Modular weights of Higgs fields are chosen to be -2. These 
choices are compatible with known results from superstring compactification where rii is -2 for 
non oscillated twisted fields of orbifold models. Anyway, this particular choice is only meant to 
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shift the degeneracy of B's and m 2 's at the unification scale in a way which is qualitatively 
consistent with our expectations from supergravity theory. We impose that masses of squarks at 
the compactification scale are all positive. This entails that: 

tan 2 #> -(n m + l) (6.11) 

where n m is the lowest modular weight. In our case tan 2 9 > 2. As a matter of fact, this is not 
strictly necessary since a negative squark mass 2 at the unification scale will be turned positive 
(mostly) by strong interaction radiative corrections, cf. (4.6). But if one assumes similar modular 
weights for the lepton sector, then (6.11) ensures that sneutrino and selectron masses are consistent 
with experimental bounds. 




Fig.l: The left mass insertion 5^ b as function of 9 (rad), 
the full line is the renormalised curve (see text). 

From (6.5), (6.6), (6.10) and (4.6), one obtains the predictions for the matrices (4.7), (4.8) in 
this particular model. It is then possible to vary 9 and plot &L,i,j as functions of 9 (see figs. 1,2 
and 3 where 9 is measured in radians). As a rule, effects of renormalisation are quite severe. On 
each plot we have depicted two curves, one without renormalisation but for the gluino contribution 
7Mq to the average squark mass 2 , compared to the renormalised one. As already discussed in the 
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introduction, the net effect of renormalisation is to reduce the value of mass insertions as the squark 
masses get a large common contribution proportional to the gluino mass [12,13,14]. Otherwise, 
the effect of renormalisation is small for the (ds) insertion whereas it is large for the other two. 
Notice that the chosen range for 9 is such that mass insertions are within the phenomenological 
bounds discussed in [8, 10, 11]. 




Fig. 2 The mass insertion <5f b as function of 6 (rad), 
the full line is the renormalised curve (see text). 



Because of the recent measurements of the b — ► sj transition, let us consider the helicity 
flipping analogue of (4.7), (4.8). They are given by the flavour mixing elements in the analytic 
sector of the mass 2 matrix. Let us define, 



~ D K(B%m D +m D B?)Kj 



and use (6.5), (6.6) and (4.6) to obtain the following numerical expressions: 



(*-L=-<- 0Scot « +0 - 9 >i 

S?S\ = -.04cot«-?^ 



Jdb M 

pf fl ) a6 = -(-06cotfl + .015)^. (6.13) 
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which are all below the experimental bounds. 




Fig. 3 The mass insertion 5% s as function of 9 (rad), 
the full line is the renormalised curve (see text). 

Since we are interested in renormalisation effects, we have concentrated on the quark sector. 
Indeed, quantum corrections are much smaller in the leptonic sector where the large QCD and top 
couplings do not contribute (quantum effects can be bigger in the large tan/3 scenario). However 
the experimental FCNC restrictions are much stronger because of the limits on \i — ► and 
similar processes. Let us generalise the mixing matrix U in (6.2) to the leptons, choose the lepton 
modular weights to be those assigned to quarks in the "same family". Then we can estimate 
the analogues of (4.8) and (6.12) in the lepton sector. By imposing the experimental bounds for 
\i — ► ej on the resulting expressions for (S^en ,(SR)en and {&LR)en, one gets a strong constraint 
on the Goldstino angle, tan 2 > 10. With this restriction on 9, the values of (<5f )ij in the figures 
lie well below their phenomenological limits. 

In summary, we have considered a simple model to estimate FCNC effects including the RGE 
running, with maximal family mixing at the unification scale. The flavour dependent supergravity 
couplings is due to the moduli fields, represented in our example by a single superfield, T. Thus 
FCNC effects are proportional to cot# ( or cot 2 9), which measures the supersymmetry break- 
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ing component along the flavour dependent direction. The resulting FCNC effects are small for 
two reasons: (i) The universal gluino contribution to the squark masses is very large, a model- 
independent feature, (ii) The contributions of the moduli sector to the scalar mass 2 are proportional 
to (rij + 1) cot 2 9 which are basically non-positive if the modular weights rii < — 1 as usually found 
in orbifold models: this fact sets a physical upper limit on the supersymmetry breaking component 
along the moduli direction. Expressed as a bound on cot# it restricts the FCNC effects to the 
relatively small values in the figures, even in the maximal mixing model discussed here. 

This negative contribution to the squark mass 2 may be considered as an artefact of the model, 
or more generally, of flavour dependence induced by the moduli couplings. Even if this is an inter- 
esting natural mechanism to generate non-universal supergravity couplings, it is certainly possible 
to contemplate other possibilities. However, unless quarks and leptons exhibit very different pat- 
terns, there is another general constraint on FCNC effects. Indeed, whatever mechanism generates 
flavour dependence, the mixing effects in the lepton sector are severely restricted by experiments. 
Hence supersymmetry breaking auxiliary fields with flavour dependent couplings are bounded, as 
in the example here above, and quark flavour changing effects are expected to be well below the 
experimental limit. But, a priori, quarks and leptons could couple to supergravity in quite different 
ways. 

7. CONCLUSION 

Within the supergravity framework, most of the unitary transformations needed to diagonalise 
the quark masses acquire a physical meaning. The RGE are then needed to make the connection 
with low energy flavour mixing phenomena. The systematic study presented here, shows how the 
quark mass hierarchies helps simplifying results. This is so even if the tree-level mixing is large. 
The effects of the running down to low energies of the squark mass matrices is sizeable in many 
respects, and the factorisation of the A-terms into Bl and Br is a useful tool. Nevertheless, some 
uncertainties remain due to the dependence on the unknown analogue of the CKM matrix for 
right-handed quarks. 

Approximate discrete symmetries could be invoked to justify that the large flavour mixing in 
the Yukawa couplings would disfavour non-universality in the Kahler potential. An elegant analysis 
[22] shows how ad hoc broken continuous symmetries could naturally control non-universality 
patterns in the squark masses and prevent FCNC problems. It would be relevant to study how 
natural these symmetries are in the context of supergravity. 
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